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D ' Abstract. It has recently been noticed that transverse momenta (pr) distributions 

observed in high energy production processes exhibit remarkably universal scaling 

behaviour. This is the case when a suitable variable replaces the usual pr- On the 
vj ' other hand, it is also widely known that transverse momentum distributions in general 

r-^... , follow a power-like Tsallis distribution, rather than an exponential Boltzmann-Gibbs 

OO ' one, with a (generally energy dependent) nonextensivity parameter q. Here we show 

J^ i that it is possible to choose a suitable variable such that all the data can be fitted 

by the same Tsallis distribution (with the same, energy independent value of the q- 

parameter). They thus exhibit q-scaling. 
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Almost fifty years ago Hagedorn developed a statistical description of momentum 
spectra observed experimentally jT]. It predicts an exponential decay of differential cross 
sections 

-$ — (-f) « 

for transverse momenta, whereas in experiments one observes non-exponential behaviour 
for large transverse momenta. Hagedorn then proposed the 'QCD inspired' empirical 
formula describing the data of the invariant cross section of hadrons as a function of pt 
over a wide range [2]: 

E^ = C.fl + ^V"^|"^t'^) f°^^^^°' (2) 

d^P \ Po) I f^y forpT^oo, 

with C, Pq and a being fit parameters. This becomes pure exponential for small px and 
pure power law for large p^lil. 

When looking at pr spectra ([T]) of secondaries produced in high energy multiparticle 
processes, it is commonly assumed that the temperature T of the hadronizing system, 
when treated as a statistical ensemble, can be connected with the observed mean 
transverse momentum (pr)- Usually, however, the system is far from thermal 
equilibrium and the individual event temperature T cannot correspond to the mean 
transverse momenta. The temperature fluctuates from event to event (or also in the 
same event). Such a situation is described by a nonextensive generalization of statistical 
mechanics proposed quite some time ago [5]. There is one new parameter, q, in addition 
to the temperature T, and the main formula of interest here is the Tsallis distribution, 

K (pt) = C, • [l - (1 - g)^] "^ ^ h (pt) = C, ■ exp (-^) , (3) 

(where Cg is a normalization constant). This coincides with Eq. (|2]) for 

1 T 

a = and po = -. (4) 

g — 1 9 — 1 

This approach has been shown to be very successful in describing very different physical 

systems [5] . Among them are also multiparticle production processes of a different kind 

(see P,[7] for recent reviews). The basic conceptual difference between ([2]) and ([3]) is in 

the underlying physical picture. In ([2]) the small pt region is governed by soft physics 

possibly described by some unknown nonperturbative theory or model, and the large 

Pt region is governed by hard physics believed to be described by perturbative QCD. In 

([3]) , the nonextensive formula is valid in the whole range of pt and it does not claim to 

originate from any particular theory. It is just a generalization of the usual statistical 

mechanics and merely offers the kind of general unifying principle, namely the existence 

of some kind of complicated equihbrium (or steady state) involving all scales of pt, 

which is described by two parameters, T and q. The temperature T describes its mean 

properties and the parameter q, known as the nonextensivity parameter, describes action 

I Actually this QCD inspired formula was proposed earlier in [3j|4]. 
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of the possible nontrivial long range effects believed to be caused by fluctuations (but 
also by some correlations or long memory effects) [5]. In fact, it was shown in |8j that 
q is directly connected to the variance of T, 

Var{T) 



q=l + 



(5) 



<T>2' 

and therefore describes its intrinsic fluctuations. For q = I one recovers the usual BG 
distribution (Eq. ([1])). In other words, the widely used thermal bath concept fails to 
satisfy conditions allowing us to introduce the notion of thermal equilibrium in the BG 
sense: it is always finite and can hardly be considered homogenous. In fact, in many 
cases it only occupies a fraction of the allowed phase space or even has a fractal-like 
structure. In such cases, a heat bath cannot be described by a single parameter T. One 
has to extend the parameter space to account for all effects mentiond above and use Eq. 
([3]) with a new additional parameter 
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Figure 1. (Color online) (a) Transverse momenta distributions considered by us. (6) 
The corresponding values of the parameter q obtained from Tsallis fits. Data are from 
compilation by Wibig [19], NA49 [17], UAl [18] and CMS [13]. 



In what follows we concentrate on data on pt distributions of secondaries produced 
in p + p and p + p collisions at selected energies (covering the wide range of available 
energies at roughly the same distance in logarithmic scale): 7 TeV and 0.9 TeV from 
CMS [13], 200 GeV from UAl [18] and 17.3 GeV from NA49 [H], cf. Fig. [Di. They 



can be fitted at each energy by using Eq. ([3]) with /i(pr) = dN/dp\. From these fits 
one finds values of q for different energies, q{s) shown in Fig. [Tjo. It can be represented 
by q{s) = 4/3 — 0.4 (yi) ' (full line). The g(s) from Fig. [T]d can be translated to 

§ Cf. [6l [7] for further references concerning specific applications of this approach to hadronic and 
nuclear physics in last decade. Recent examples of power-like distributions fitted by nonextensive 
Tsallis formula are provided in [9] [10] , by PHENIX [U, and STAR J^J experiments at RHIC and 
CMS aZ, ATLAS [H] and ALICE [E] experiments at LHC. Finally, the possible QCD origin of 
such fluctuations and/or correlations could probably be traced back to the nonperturbative QCD, for 
example, |16j . 
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Figure 2. (Color online) (a) Experimental values of mean transverse momenta, 
(pt(s)), of charged particles produced in p + p and p + p collisions [13l [171 HI] (cf- 
Fig. [T^ for corresponding distributions at selected energies). (5) The resulting T{s) as 
given by Eq. ([6]). 



energy dependence of the temperature, T(s), with the help of {pt{s)), 

T 

(Pt) 



4- 3,' (^' 

Using {pt{s)) as evaluated experimentally, cf. Fig. [2^ (it can be parameterized by 
{pris)) = 0.235 (y/s) ' ) one gets the T(s) presented in Fig. [2b, which can be 
represented by T{s) = 0.141 {^/s) ' . To summarize: all pr data considered here 
can be fitted with a Tsallis formula, Eq. ([3]), by using an energy dependent parameter 
q{s) and T(s). 

Recently there have been attempts to use an energy independent q to describe 
dN/dpT and to check for possible scaling behaviour in pt [20]. It was found that to this 
end one has to use a p^-dependent form of the nonextensivity parameter, 

go- (go- l)^(Pr) 



(l{Pi 



(7) 



l-{q,-l)e{pT) 

where go = 1.12 ± 0.06 and 9 {pr) = 21og [log (1 + k^t)] with k = 0.013 GeV^\ It 
was shown that any experimentally accessible g depends on pt growing from g = 1 for 
Pt = and attaining g = go for pt = 9/k. 

The distributions dN/dpx shown in Fig. [TK differ for different energies. However, 
as shown in |211 [221 123], one can find a single scaling function F{t), independent 
of energy, and a suitable scaling variable r such that one observes a scaling: 
h{pT, y/s) -^ F {t = f {pt, y/s)), analogous to Feynman or KNO scaling [21]. Prompted 
by geometrical scaling behaviour found in deep inelastic scattering data, the following 
universal variable has been proposed 
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with Ql^t {pt) = Ql 



W 



(8) 
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Figure 3. (Color online) (a) Data for transverse momentum distributions for different 
energies [T71 [TSJ [T3] plotted by using the scaling variable p'j, defined in [23^ and fitted 
using Eq. ([3|) . (6) The ratio data/fit for the results presented in (a) . 



In this variable all data lie on a single curve, 

dN 1 



Q' 



Ft) 



(9) 



dpT ^0 

with F{t) being some energy independent universal function (cf. [23] for details) and A 
is a parameter. Actually, to get good agreement with all available data, A has to depend 
on pt- The best fit, see Fig. El is obtained with A = XeffiQ) = 0.13 + 0.1 {Q'^/10f^\ 
where Q = 2pT (cf., Eq. (11) of [23]). The variable p'rp used in Fig. [3] was obtained 
by demanding that pt at energy W should be connected with p'j, at energy W via the 
following relation (cf. [23] for details). 



Pt = Pt 



'w 



X 

X+2 



(10) 



Fig. |3] shows the corresponding results together with a Tsallis fit performed using Eq. 
(^ with C = 14.0, q = 1.121, T = 0.18 GeV. 

Here we would like to check whether already considered data show a kind of q- 
scaling as seen from the perspective of Tsallis statistics (and, if so, in what variable). 
In other words: is it possible to find, in the framework of the nonextensive statistics, 
a variable (different from the p'j, above) which would scale the dN/dpT distributions? 
And, is the parameter q = 1.121 from the Tsallis fit in Fig. [3|i already universal, or else 
can one also have such scaling behaviour for some other value of the parameter q using 
a different scaling variable? 

Prompted by KNO scaling as observed in multiplicity distributions [23] we first plot 
data from Fig. [1^ using the scaled transverse momentum variable, 
Pt 



{PtY 



ir 



As seen in Fig. HJ already this variable seems to be nearly satisfactory, except for the 
largest LHC energies. Agreement with data can be further improved by using Tsallis 
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Figure 4. (Color online) Data for transverse momentum distributions for different 
energies [ITl [181 US] plotted by using a scaling variable z as defined in Eq. (|lip . 



distribution, hq{u) 



u 



as given by Eq. 
Pt 



), in which Pt/T is replaced by m/mq, wher 



{Pt) -b-pT 

Using {pt{s)) taken from an experiment as shown in Fig. [2^, and an energy dependent 
coefficient b{s) = —0.0397 + 0.08 (a/s) ' , one can reasonably well fit the data, cf.. 
Fig. [5l with a constant, energy independent value of g = 1.172 (and with C = 79.4 and 
Mo = 0.17). Notice that our result is essentially of the same quality as that obtained 
from the geometrical scaling prescription proposed in [23] (cf.. Fig. [3]). Both in Fig. 
|3] and in Fig. [5] there are deviations from Tsallis distributions. Essential here are the 
differences between different energies. To show this in both cases we evaluate the ratios 



d^N 



2'KpTdpTdy 



(13) 
for different energies y^, which 



^ 7 (pt, v^ = 7 TeV) 

of experimental distributions, / {pr, a/s) 
are shown in Fig. |6l 

To justify using the variable u as defined in Eq. (IT2l) . note that one can write 
u/uq = pr/Teff where T^ff is an effective temperature, 

Teff = To + T^ (pt) ; with Tq = uq ■ (pr); T^ = -b ■ Uq ■ pr- (14) 

This temperature could be related to the possible pt transfer, additional to that resulting 
from a hard collision, perhaps proceeding by a kind of multiple scattering process. 



Notice that u is just a power series in the scaling variable z defined in Eq. ([TT|) . 
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Figure 5. (Color online) (a) Data for transverse momentum distributions for diflFerent 
energies [T71 [TSJ [T3] plotted by using the scaling variable u defined by Eq. ([T^ and 
fitted by Eq. ([3|) in the variable u. (6) The ratio data/fit for results presented in (a). 
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Figure 6. (Color online) Ratios of transverse momenta distributions at energies 
^/l = 0.9 and 0.2 TeV with reference to the distribution at ^/s = 7 TeV, expressed in 
variable p'rp at (a) and variable u at (b). 



similar, in a sense, to that proposed on a different occasion in [25]. It is therefore 
not necessarily connected with thermodynamics. In fact such Tg// also occurs in a 
description of the growth of the so called complex free networks if one associates pt with 
the number of links [26]. When looking at hadron production from the perspective of 
stochastic networks [27j one can argue that the power law seen in transverse momenta 
spectra means that hadronization can be viewed as a process of formation of some 
specific network taking place in the environment of gluons and quark-antiquark pairs 
formed during the hadronization process. In this case their actual original energy- 
momentum distributions would be of secondary importance in comparison to the fact 
that, because of their mutual interactions, they connect to each other and that this 
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process of connection has its distinctive dynamical consequenceqj]]. 
More formally, following |26] observe that, whereas 



dfjx) 
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~fi^) 



fix) 



^exp 
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the X- dependent T in the form, 

T ^ T{x) = To + (g - l)x 
results in a Tsallis distribution: 

df{x) 1 



(15) 



(16) 



dx 



To + (g - l)x 



fix) 



fix) 
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Tn 
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Q)7f- 
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.(17) 



When fitting data with this distribution, one encounters the necessity to use in Eq. fITTj) 
the s-dependent q and Tq: 

q = q(s) = go + q'{s) and Tq = To(s) = T^ + T'(s). 

To compensate for this s-dependence, one can modify the temperature in Eq. (p^ 
example by allowing for an x-dependence: 

To -^Teff{x) =TQ-bx, 
where the parameter b can be s-dependent. Returning to Eq. (TT6l) . one now has 

T{x) = Tg + (g — l)x — bx 
and that, solving the present form of Eq. (IT71) . one finds 



(18) 
, for 

(19) 
(20) 



^0 


Identifying now: xb{s) = T'{s) 
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(21) 



and b{s) = g'(s), one obtains an energy independent 
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(22) 



In reality the situation is more complicated because here the x-dependence was 
introduced to T on the level of the distribution function, not in the differential equation. 

^ The possible line of reasoning is as follows: Suppose we start with some initial state consisting of a 
number no of already existing {qq) pairs (identified with vertices in the network). We add to them, in 
each consecutive time step, another vertex (a new {qq) pair), which can have fcg possible connections 
(links in the network language) to the old state. Assume that quarks are dressed by interaction with 
surrounding gluons and therefore "excited" and that each quark interacts with k other quarks (has 
k links). Assuming further that the "excitation" of a quark is proportional to the number of links k 
(which is proportional to the number of gluons participating in "excitation" , i.e., existing in the vicinity 
of a given quark) , the chances to interact with a given quark grow with the number of links k attached 
to it. The new links will be preferentially attached to quarks already having large k. This corresponds 
to building up a so called preferential network, which evolves due to the occurrence of new {qq) pairs 
from decaying gluons. 
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Nevertheless, it seems that with such manipulations one can expect a distribution of 

the form 

1 



f{x) oc 



1 - (1 - g, 
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^-leff 



(23) 



in which the s-dependence should be noticeably reduced. This leads us to the variable 
u introduced in Eq. flT2|) . 
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Figure 7. (Color online) (a) Data for transverse momentum distributions for different 
energies ^7\ [HI [T^ plotted by using the scaling variable u defined by Eq. ([T2|) with 
& < 0. (6) Ratios of transverse momenta distributions at energies y/s — 0.9, 0.2 and 
0.0173 TeV with reference to the distribution at y/s — 7 TeV . 
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Figure 8. (Color online) (a) Data for transverse momentum distributions for central 
A+A collisions at different energies. Data for Pb+Pb at 2.76 TeV comes from ALICE 
[H] and CMS HH], data for Au+Au at 0.2 TeV come from PHENIX [55] and for Pb+Pb 
at 0.0173 TeV from NA49 [30]. (b) The same plotted by using the scaling variable u 
defined by Eq. ([12]) with & < 0. 



Notice that, whereas for the choice of b{s) used above one has 6 > and Tg// was 
decreasing with px (increasing, in the network approach discussed above, action of the 
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preferential attachement) one can also choose a parametrization for which 6 < and 
Teff increases with p^ cancelling effects of "preferential attachment". For example, for 
b{s) = —0.109 + 0.115 (y^) ' , one obtains a scaling of distributions in the variable 
u with the quasi-exponential form of the function h{u) (i.e., in a Tsallis distribution 
with q close to 1) . The corresponding results are shown in Fig. [TJ together with the 
corresponding Tsallis fit for q = 0.955 (and C = 5, Uq = 0.6). 

Actually, the same kind of scaling is also possible for A + A collisions for b{s) = 
—0.052 — 0.0002 (y/s) ' . In Fig. [8] one can see distributions oi px for central collisions 
(0 — 5% centrality) for different data together with the Tsallis fit for q = 0.9999 (and 
C = 2800, uo = 0.32). 

To summarize: it is possible to fit all available data on pt distributions using some 
universal, energy independent, parameter q. Therefore there is a possibility of g-scaling. 
This can be done by choosing a variable u defined in Eq. ( IT2l) in the distribution h{u) 
given by Eq. ([3]). Scaling can be achieved either by increasing q to q = 1.172 (for 6 > 0, 
cf.. Fig. [5]) or by decreasing it to g ~ 1 (for 6 < 0, cf.. Fig. O in this case Eq. ([3]) almost 
coincides with Eq. ([I]))lJ. The observed Tsallis distributions do not necessarily indicate 
thermalization of the system considered (there are numerous examples of non-thermal 
sources of Tsallis distributions, cf. [7]). The possible explanation we propose is based on 
the description of a hadronization process in analogy with complex free networks |27j . 
Alternativelly, one can interpret Eq. (ll7p as a stationary solution of the Fokker-Planck 
equation corresponding to a Langevin equation with multiplicative noise with nonzero 
mean value {r](t)) = 1 — q [31] (cf. also [7]). The possible connection with QCD based 
ideas [25] is also indicated (but this would demand special attention, which is outside 
of the scope of this paper). 
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